Magnetic monopoles and synthetic spin-orbit coupling in Rydberg macrodimers 
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We show that sizeable Abelian and non-Abelian gauge fields arise in the relative motion of two 
dipole-dipole interacting Rydberg atoms. Our system exhibits two magnetic monopoles for adiabatic 
motion in one internal two-atom state. These monopoles occur at a characteristic distance between 
the atoms that is of the order of one micron. The deflection of the relative motion due to the 
Lorentz force gives rise to a clear signature of the effective magnetic field. In addition, we consider 
non-adiabatic transitions between two near- degenerate internal states and show that the associated 
gauge fields are non-Abelian. We present quantum mechanical calculations of this synthetic spin- 
orbit coupling and show that it realizes a velocity-dependent beamsplitter. 
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Gauge theories represent a cornerstone of modern 
physics and play a prominent role in classical and quan- 
tum electrodynamics, the standard model of elementary 
particle physics and condensed matter physics. In view 
of the importance of this concept tremendous effort has 
been made to create artificial gauge fields for neutral 
atoms |lMl2| and to investigate the resulting atom dy- 
namics in the quantum regime. In these schemes, engi- 
neered light-matter interactions cause neutral atoms to 
behave like charged particles in an electromagnetic field. 

Artificial gauge fields allow the simulation of theoreti- 
cal models that are otherwise inaccessible. For example, 
the realization of magnetic monopoles affecting the rela- 
tive nuclear motion of diatomic molecules has been dis- 
cussed in [13|. However, gauge field effects in molecules 
are usually very small since they arise from terms that 
are neglected in the Born-Oppenheimer approximation 
(BOA) [14| which is very well satisfied in many molecu- 
lar systems. In addition, the experimental observation of 
these effects is considerably hampered by the small size 
of conventional molecules. 

Recently, extremely large molecules comprised of two 
Rydberg atoms with non-overlapping electron clouds 
have been proposed [15l4l7| and observed [18J. Typical 
internuclear spacings exceed 1 /im, and thus the experi- 
mental observation of Rydberg- Rydberg correlations be- 
comes feasible 19I422J. These so-called macrodimers in- 
teract via well-understood and controllable dipole-dipole 
potentials. Importantly, the validity of the BOA cannot 
be established via the mass ratio of nuclei and electrons 
for these systems. 

Here we show that dipole-dipole interacting Rydberg 
atoms can exhibit Abelian and non-Abelian gauge fields 
that influence the quantum dynamics of the relative 
atomic motion substantially. In contrast to the Rydberg 
macrodimer proposal in [17j, the system in Fig. [I] is dis- 
tinguished by an asymmetric Stark shift of the Zeeman 
sublevels. We find that this broken symmetry gives rise 
to magnetic monopoles if the system evolves adiabatically 
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FIG. 1. (Color online) (a) System geometry of two dipole- 
dipole interacting Rydberg atoms. R is the relative position 
of atom 2 with respect to atom 1, and p is the distance of atom 
2 from the z axis. The z direction is distinguished by a DC 
electric field E. (b) Internal level structure of each Rydberg 
atom. The Stark shifts 5 and A are negative. We assume 
S ^ A. The dipole transitions indicated by blue dotted, black 
solid and red dashed lines couple to a~ , tt and o~ + polarized 
fields, respectively. 



in one internal two-atom state. These monopoles occur 
at a characteristic distance between the atoms that is of 
the order of one micron. The Lorentz force associated 
with the magnetic field near a monopole results in a size- 
able deflection of the relative atomic motion. This effect 
can be interpreted in terms of an exchange between or- 
bital and internal spin angular momentum as the internal 
molecular state changes while the atoms move. Moreover, 
we investigate non-adiabatic transitions between two in- 
ternal two-atom states and find that the associated gauge 
fields are are non-Abelian. This synthetic spin-orbit cou- 
pling creates a coherent superposition of two spatial con- 
figurations of the atoms. We expect that our findings are 
relevant for other dipole-dipole interacting systems like 
polar molecules and magnetic atoms. 

The geometry of the two atom system under consid- 
eration is shown in Fig. [Ha). In order to account for 
the azimuthal symmetry of the system, we express the 
relative position R of atom 2 with respect to atom 1 in 
terms of cylindrical coordinates, R = (pcos(/>, psin^, z). 
We omit the center-of-mass motion which is uniform and 



investigate the relative motion of the two dipole-dipole 
interacting Rydberg atoms. The Hamiltonian of this sys- 
tem is given by 



H = H R + H s + V dd , 



(1) 



where #r = p 2 /(2fi) is the kinetic energy of the rela- 
tive motion and fi is the reduced mass. H$ describes 
the internal levels of the two uncoupled atoms and Vdd 
is the dipole-dipole interaction [17J. In each Rydberg 
atom we consider two angular momentum multiplets as 
shown in Fig.[TJb). The lower nsi/ 2 states have total an- 
gular momentum J = 1/2, and the excited multiplet is 
comprised of np 3 / 2 states with total angular momentum 
J = 3/2. We specify the individual atomic states \£rrij} 
by their orbital angular momentum £ and azimuthal total 
angular momentum rrij. A DC electric field E in the z 
direction defines the quantization axis and gives rise to 
Stark shifts of the magnetic sublevels. We assume that 
the Stark shifts are different in the rrij > and rrij < 
manifolds, which could be achieved, e.g., by inducing ad- 
ditional AC stark shifts. For simplicity we focus on the 
level scheme shown in Fig. H^b), where the asymmetry is 
characterized by the ratio A/5 of the Stark shifts 5 and 
A. The relevant subspace of two- atom states is spanned 
by the N = 16 nsnp states where one atom is in a nsi / 2 
state and the other in a np 3 / 2 state. For every value of 
R we introduce a set of orthonormal eigenstates of the 
Hamiltonian H$ + Fdd, 



(H s + v dd )\MR)) = ei(R)\MR)), 



(2) 



where Ci(R) are the corresponding eigenval- 
ues. With these definitions, the full quantum 
state of the two-atom system can be written as 
l*> = Eti/d 3 i? ai(R)\^i(R)) \R). Next we 
assume that the dynamics is confined to q eigenstates of 
Hs + Vdd? i- e -? there may be non-adiabatic transitions 
within the first q eigenstates, but transitions to other 
states \ipi(R)) (I > q) can be neglected. We follow the 
procedure described in pj-y, |23| and derive from Eq. (pp) 
an effective Schrodinger equation for the wavefunctions 
ol = (ai,...,a q ), q< N, 



ihdtcx 
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(p - A) 2 + V + $ 
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(3) 



Equation (J3]) is equivalent to the Schrodinger equation 
of a charged particle in an electromagnetic field, charac- 
terized by the vector potential A and scalar potential $. 
Here V and A are q x q matrices whose matrix elements 
for k, I < q are given by 

V H = 5 kl e k (R), A kl = ih(MR)MMR)), (4) 

and Ski is the Kronecker delta. We find that the impact 
of the scalar potential $ on the presented results is neg- 
ligible, and thus omit it in the following. The Cartesian 
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FIG. 2. (Color online) (a) Potential curve e\ in the x — y 
plane as a function of scaled internuclear spacing R/Ro . The 
parameters are A = — 1.5\5\ (dashed black line) and A = 
—3\5\ (red solid line), (b) Potential curve e± in the x — y 
plane for A = — 3\S\. 



components B^ (i G {1,2,3}) of the artificial magnetic 
field are defined as 

B (i) = \e tkl F^ kl \ (5) 

F (ki) = dkA (i) _ diA {k) _ | j A ( fc ) > A (i) j ^ (6) 

where Siki is the Levi-Civita tensor and we employed Ein- 
stein's sum convention. The q x q matrices A^ describe 
non-Abelian gauge fields if the commutator [A^,AO] 
is different from zero. The magnetic field gives rise to 
a Lorentz force which is proportional to the velocity 
v = (p — A) / [i of the relative motion. 

The eigenstates |^) and eigenvalues q of the Hamil- 
tonian Hs + Kid m Eq. ([2]) can be obtained numerically. 
Here we focus on one particular potential curve that ex- 
hibits a potential well such that the two dipole-dipole 
interacting Rydberg atoms can form a giant molecule. 
This potential curve is labelled by e\ and shown in 
Fig. [2ja) for different ratios of A/5. The potential min- 
imum occurs roughly at the characteristic length Ro = 
(|P| 2 / (4^60^1(5 1) 1 / 3 denoting the distance where the mag- 
nitude of the dipole-dipole interaction equals the Stark 
splitting h\5\ [17[, and \V\ is the reduced dipole ma- 
trix element of the ns «-»• np transition [17| . Due to the 
azimuthal symmetry of the system a donut shaped po- 
tential well arises in the x — y plane [see Fig. [2fb)]. The 
dependence of e\ on x and z is displayed in Fig. [3fa), 
showing that the width of the potential well in z direction 
is comparable to its width in the x — y plane. Note that 
a potential well with similar features but for 5 = A < 
was reported in [17| . 

Next we consider the case q = 1 in Eq. ([3]) and consider 
the adiabatic motion in the eigenstate |^i) corresponding 
to ei. We choose the phase of |^i) such that the vector 
potential A\ = A\\ obeys the Coulomb gauge (divAi = 
0) 0. It follows that 



A™(p,z) = A 1 .e 4 



1 
P 



(MR)\Jz\MR)) (7) 



is the only non-zero component of the vector potential. 

In this equation, J z = J, 

ponent of the total angular momentum operator of the 
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Jz and Jz is the z com- 
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FIG. 3. (Color online) (a)-(c) Abelian case. We consider 
adiabatic motion in state |-0i). (a) Potential curve ei in the 
x — z plane for A = — 3|£|. (b) Artificial magnetic field B 
in the p — z plane. R± indicate the positions of the mag- 
netic monopoles. (c) Color code for the absolute value of 
the effective magnetic field in (b). (d)-(f) Non- Abelian case. 
We consider non- adiabatic motion in states \ipi) and ^2) for 
A = — 1.16 1^| . (d) Potential curves e\ and 62 in the x — y 
plane. The initial position of the Gaussian wavepacket in 
the upper well £2 is indicated by a black arrow, (e) and (f) 
show the probability density of the evolved wave functions a± 
and cx<2 at t\S\ = 220. The population in the upper (lower) 
well state |^ 2 ) (|^i)) is 52.7% (47.3%). In (e) and (f), we 
consider 39 K atoms with principal quantum number n = 20, 
\S\ = 2tt x 24.7 MHz and R = 1.28/xm. 



internal states of atom a and e^ is the unit vector in <p 
direction. The spatial variation of the vector potential 
determines the magnetic field B according to Eq. ([5]). 
We find that B is only different from zero if S 7^ A 
such that the symmetry of the system is broken. In 
contrast, for a symmetric level scheme the expectation 
value of J z in Eq. (|7|) yields zero. Figure [3fb) shows 
the magnetic field in the p — z plane for A = — 3\S\. 
The most remarkable features of B are a source and 
a drain of magnetic flux near R+ « +0.98 Roe z and 
R- « —0.98Roe z , respectively. We integrate the mag- 
netic field over a sphere S± centered at R± and evaluate 
the Chern number C± = J s± B • dS/(27r). The re- 



sult is C± = ±1, demonstrating that our system exhibits 
two magnetic monopoles on the z axis. We emphasize 
that the monopoles arise at an atomic separation roughly 
given by R§. This parameter can be controlled by the 
magnitude of the dipole-dipole interaction and the Stark 
splitting, and is typically of the order of one micron. 

Next we show that the artificial magnetic field gives 
rise to a sizeable deflection of the relative atomic motion. 
To this end, we suppose that the relative position of the 
two atoms is initially given by R\ = Ro(—1.5e z +O.Obe x ). 
The initial velocity of the relative motion is (v) = v$e z 
with ^o > such that the atoms move towards the re- 
gion of strong magnetic fields near R_. Here we treat 
the relative atomic motion classically, which is justified 
if the wavepacket associated with the relative motion is 
very well localized. Since R\ contains a small offset in the 
positive x direction, the potential curve e\ will result in a 
deflection in the positive x direction [see Fig. 03(a)] . Note 
that our choice of initial conditions and the azimuthal 
symmetry of the system imply that the motion remains in 
the x — z plane if there were no magnetic fields. However, 
the magnetic field pointing towards R- yields to a deflec- 
tion y_ in the negative y direction via the Lorentz force. 
On the contrary, the magnetic field will have the oppo- 
site effect if we mirror our initial conditions at the y — z 
plane, i.e., for Ri = Ro(— 1.5e z — 0.05^). In this case, 
the relative motion remains in the x < half plane and 
the Lorentz force results in a deflection y+ in the positive 
y direction. It follows that the difference Ay = y+ — y_ is 
a direct measure of the effective magnetic field. In addi- 
tion, it reflects the broken chiral symmetry arising from 
the asymmetric level scheme in Fig. [2(b). 

In order to obtain a quantitative description of this ef- 
fect, we derive a set of coupled equations for the mean val- 
ues (R) and (v) from the Hamiltonian H in Eq. (pQ) [261 ]. 
We consider 23 Na atoms with principal quantum num- 
ber n = 15 and \S\ = 2tt x 39.0 MHz 0. This yields 
Ro = 0.75/im, and hence the distance between the atoms 
is initially given by \Rj\ = \Rj\ « 1.13/im. Furthermore, 
we set ^o = 195ft/ (ifo/i). We neglect effects due to the 
finite lifetime of the molecule (T\S\ « 540) and thus re- 
strict the analysis to times t < T. From semi-classical 
simulations for the two initial positions Ri and Ri we 
find Ay = O.LRo for t\S\ = 238, and Ay = 0.3R at 
t\5\ = 483. It follows that the magnetic field results in a 
substantial deflection of the relative atomic motion. Note 
that our simulations allow us to confirm that the motion 
remains adiabatic at all times. 

Next we show that the vector potential A can give rise 
to a coupling between the relative atomic motion and 
internal electronic states. In order to demonstrate this 
synthetic spin-orbit coupling, we consider an additional 
potential curve 62 with corresponding state l^)- Here we 
focus on the two-dimensional setting where the motion 
is confined to the x — y plane. The two potential curves 
ei and 62 become near-degenerate for R « L33.Ro and 




FIG. 4. (Color online) (a) Potential curves e\ (solid blue line) 
and €2 (dashed blue line) in the x — y plane, (b) Imaginary 
part of A$ = [A^]* for = 0. (c) Real parts of A { $ (solid 
red line), A 22 (dashed black line) and A\ 2 = [A 21 ]* (dotted 
blue line) for = 0. (d) Matrix elements of the commutator 
C = z[A (1 \A (2) ]/h. The red solid line shows Cn = -C 22 , 
and the black dashed line represents Cn = Ci\. In (a)-(d), 
we set A = — 1.16|£|. All components of A that are not shown 
in (b) and (c) are zero. 



A = — 1.16 1 J | and are shown in Fig. HJa). In order to in- 
vestigate the quantum dynamics in the two cylindrically 
symmetric potential wells, we evaluate the vector poten- 
tial A in Eq. (|4]) by numerical means [24j . Note that A is 
now represented by a 2 x 2 matrix, where each component 
Aki is a 3-column vector. We find that the component 
A^ is zero, and the non-zero parts of A^ and A^ are 
shown in Figs.|U^b) and (c), respectively. All components 
of A are evaluated for <\> = such that A^ (A^) can be 
identified with the radial (azimuthal) component of A. 
Near the avoided crossing the off-diagonal element A\2 
can induce non-adiabatic transitions between the states 
|^i) and |^2) • The coupling strength depends on the en- 
ergy difference | e\ — e<i\ and the velocity of the relative 
motion. For a quantitative description of this synthetic 
spin-orbit coupling, we assume that the system is ini- 
tially at rest and prepared in the upper well state \1Jj2} 
[see Fig. 13(d)]. We model the wavepacket correspond- 
ing to the relative atomic motion by a Gaussian with a 
full width at half maximum of a « 75 nm centered at 
R = 1.5Ro and solve Eq. (J3J) for q = 2 in a box with ra- 
dius 2.2Ro. As the system evolves, it will oscillate in the 
upper well, and near the avoided crossing some popula- 
tion will be coherently transferred to the lower well state 
|^i). The probability densities in the two states after 
the avoided crossing has been traversed once is shown in 
Figs. [3(e) and (f). For the chosen parameters an almost 
equal superposition of the two internal states is created. 
Note that the two wavepackets experience different po- 
tentials [see Fig. Ufa)] and hence they will separate in 
space for longer evolution times. 

We emphasize that the gauge fields A^ and A^ are 



strongly non-Abelian as shown in Fig. H(d). The com- 
mutator C = i[A^ l \ A^] / h is of the same order of mag- 
nitude as the first term in Eq. (|6]), and thus the non- 
Abelian signature is significant whenever the magnetic 
field gives rise to sizeable effects in the quantum dynamics 
of the system. This opens up the possibility to study the 
rich physics resulting from non-Abelian gauge fields [28| , 
which is subject to further investigation. 

In summary, we have shown that the dipole-dipole in- 
teraction between Rydberg atoms can induce Abelian 
and non-Abelian artificial gauge fields that influence the 
relative atomic motion significantly. The experimental 
realization of our scheme could be achieved in optical 
lattices where the lattice constant matches the desired 
initial separation of the atoms. Alternati vely , one could 
start with a similar setup as described in [2l|, |22j , where 
the dipole-dipole interaction between two individual Ry- 
dberg atoms trapped in optical tweezers was investigated. 
The optical potentials allow one to control the initial po- 
sition of the atoms before they are excited to the diatomic 
nsns state via laser fields. A subsequent microwave field 
prepares the system in the desired nsnp state |^i) or 
1^2) • I n addition, the optical trapping potentials could 
transfer linear momentum to the atoms before the ex- 
citation to the Rydberg states occurs. Our calculations 
for the deflection in the monopole field were carried out 
at zero temperature. By considering a thermal velocity 
distribution, we estimate that the deflection pattern will 
be washed out if the temperature exceeds approximately 
100 nK. These temperatures are routinely achieved in op- 
tical lattices and dipole traps [29|. Finally, the observa- 
tion of the relative atomic motion requires measurements 
of the density-density correlations of the two Rydberg 
atoms. Such measurements have been performed by ion- 
ization of the Rydberg atoms [19[ and by de-excitation 
to the ground state followed by advanced imaging tech- 
niques [20[. We thus believe that the experimental ob- 
servation of the deflection in the monopole field and the 
splitting of the motional wavepacket is feasible with cur- 
rent or next-generation imaging techniques. 
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